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Q I To study the peculiarities of the Galactic spiral density wave, we have analyzed the 

Q ■ space velocities of Galactic Cepheids with proper motions from the Hipparcos catalog 

■ and line-of-sight velocities from various sources. First, based on the entire sample of 
185 stars and taking Rq = 8 kpc, we have found the components of the peculiar solar 
velocity {uq,Vq,Wq) = (7.6,11.6,6.1) ± (0.8,1.1,0.6) km s~^, the angular velocity of 

I Galactic rotation fio = —27.4 ± 0.6 km s^^ kpc~^ and its derivatives f2g = +4.07 ± 0.21, 

^ ■ km s~^ kpc~^ and i7g = —0.83 ± 0.17, km s~^ kpc~^, the amplitudes of the velocity 

■ perturbations in the spiral density wave /r = —6.7 ± 0.7 and fg = 3.5 ± 0.5 km s~\ the 
Qh, pitch angle of a two-armed spiral pattern {m = 2) i = —4.5 ± 0.1° (which corresponds 
O ! to a wavelength A = 2.0 ± 0.1 kpc), and the phase of the Sun in the spiral density wave 

Xo = ^191 i 5°. The phase Xe has been found to change noticeably with the mean age 
ci ' of the sample. Having analyzed these phase shifts, we have determined the mean value 

of the angular velocity difference Qp — Q, which depends significantly on the calibrations 
used to estimate the individual ages of Cepheids. When estimating the ages of Cepheids 
^ ! based on Efremov's calibration, we have found — fiol = 9 ± 2 km s~^ kpc~^. The ratio 

I of the radial component of the gravitational force produced by the spiral arms to the total 

gravitational force of the Galaxy has been estimated to be fro = 0.04. 
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INTRODUCTION 



X. 

\-i . Data on various objects are used to determine the Galactic rotation and spiral struc- 

- - -' ture parameters. Classical Cepheids, whose distances are determined from the pe- 

riod. luminosity relation, are the most important objects for solving this problem (Feast 
and Whitelock 1997; Mishurov and Zenina 1999; Mel'nik et al. 1999; Rastorguev et al. 
1999). 

Cepheids are distributed in a fairly wide solar neighborhood {t^5 kpc). At present, 
highly accurate data are available for quite a few of them (about 200) to determine their 
space velocities. 

An important parameter of the Galactic spiral structure is the spiral pattern speed 
Qp. Mishurov et al. (1979) proposed a method for estimating this quantity using the 
correction factors -FJ^-*(x) and -Fj^-'(x) and found ilp = 19.1 ± 3.6 km kpc~^ from 
Cepheids. Popova (2006) obtained an estimate of Qp = 21.7 ±2.8 km s^^ kpc^^ also from 
Cepheids but by a different method based on analysis of the "Ini? — 6'" plane. At present, 
however, we have no firm confidence in the accuracy of Qp. For example, according to the 
review by Gerhard (2011), the values of Qp obtained by different authors lie within the 
range from 15 to 30 km s~^ kpc~^. 



The fact that our analysis of the radial velocities (V^) for Galactic masers (Bobylev 
and Bajkova 2010; Stepanishchev and Bobylev 2011) and 0B3 stars (Bobylev and Bajkova 
2011) showed a significant difference between these two samples of young objects at the 
Sun's phase Xq of about 50° served as one of the incentives to perform this work. This 
suggests that this difference is due to the observed difference {Qp — Q) ■ t. If accurate 
estimates of the individual or group ages for stars were available, then information about 
Qp could be directly extracted from the observed difference of the Xq values. Cepheids 
are such stars with well-known age estimates. 

The goal of this paper is to determine the Galactic rotation parameters and Galactic 
spiral density wave parameters from the space velocities of Cepheids. To estimate the 
spiral pattern speed, we suggest using a direct method based on analysis of the change in 
the Sun's phase with time. For this purpose, wc produce and investigate three samples 
of Cepheids with different mean ages. 

1 DATA 

We used data on ^240 classical Cepheids with proper motions mainly from the Hipparcos 
catalog (van Leeuwen 2007) and line-of-sight velocities from various sources. The data 
from Mishurov et al. (1997) and Gontcharov (2006) as well as from the SIMBAD and 
DDO databases served as the main sources of line-of-sight velocities for the Cepheids. 
For several long-period Cepheids, we used their proper motions from the TRC (Hog et al. 
2000) and UCAC3 (Zacharias et al. 2009) catalogs. 

To calculate the Cepheid distances, we use the calibration from Fouqu et al. (2007). 
The calibration of the period. luminosity relation for Cepheids from Berdnikov et al. (2000) 
is also well known. Only nine Galactic Cepheids were used for the calibration of Berdnikov 
et al. (2000). However, all of them are members of open clusters and pulsate in the 
same mode. The distance calibration for these open star clusters was first thoroughly 
studied. Fouqu et al. (2007) refined the period. luminosity relation for Cepheids using 59 
calibration stars with their parallaxes measured by various methods. According to Fouqu 
et al. (2007), My = -1.275- 2.678 -IgP,, where the period is in days.l Given My,, taking 
the period-averaged apparent magnitudes iVi and extinction Ay = 3.23- E{B — V) mainly 
from Acharova et al. (2012) and, for several stars, from Feast and Whitelock (1997), we 
determine the distance r from the relation 

^^-^^Q-0.2{Mv-V-b + Av)^ (1) 

and then assume the relative error in the Cepheid distances determined by this method 
to be 10%. The hst of Cepheids (their numbers according to the Hipparcos catalog) from 
Feast and Whitelock (1997) included in our sample is given in Table 1. The distances 
from Fouqu et al. (2007) and Berdnikov et al. (2000) are compared in Fig. 1. We see 
from this figure that the differences between them are insignificant, although the scale 
from Berdnikov et al. (2000) is slightly shorter. For checking, we determined all of the 
sought-for parameters using the distances r from Berdnikov et al. (2000) but found no 
significant differences in the parameters being determined. 

We divided the entire sample into three parts, depending on the pulsation period, 
which rcfiects well the mean Cepheid age. Several calibrations proposed to estimate the 
mean Cepheid age are known. We use two of them: first, the theoretical calibration from 



Table 1: Cepheids from the list by Feast and Whitelock (1997) included in our sample 



HIP 


HIP 


5138 


78797 


5846 


78978 


34527 


87072 


39144 


89013 


42492 


91366 


42926 


91613 


50615 


91738 


51653 


93399 


52380 


98376 


53867 


104002 


56991 


106754 


57649 


116556 



Bono et al. (2005), 

log ti =8.31 -0.67 log P, (2) 

for the fundamental period of Cepheids with a mean metallicity of 0.02 typical of Galactic 
stars; and, second, the calibration from Efremov (2003), 

logi2 = 8.50 - 0.65 log P, (3) 

obtained by analyzing Cepheids in the Large Magellanic Cloud. 

We rejected the double Cepheids and short-period Cepheids classified as DCEPs, which 
pulsate in the first overtone and, therefore, their distances have a low accuracy. The main 
errors in the space velocities of Cepheids are associated with the errors in their proper 
motions. This is especially clearly seen in the velocity components W. Therefore, we used 
a constraint on the absolute value of the residual (after the subtraction of the Galactic 
rotation parameters) velocity, |V{7yvt^| < 50 km s~^. All of the remaining Cepheids are no 
father than 5 kpc from the Sun. 



2 THE METHODS 
2.1 Simultaneous Solution 

This method for determining the kinematic parameters consists in minimizing a quadratic 
functional F: 

min F = E^^^«^(V7-V7)2 

provided the fulfilment of the following constraints derived from Bottlinger's formulas 
with an expansion of the angular velocity of Galactic rotation Q into a series to terms 
of the second order of smallness with respect to v/Rq, and with allowance made for the 
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Figure 1: Cepheid distances found from the calibrations of Berdnikov et al. (2000) and 
Fouqu et al. (2007). 



influence of the spiral density wave: 



Vr — —Uq cos h COS I 

— Vq cos fe sin Z — ^0 sin b 
+ Ro{R — Rq) sin I cos bQ^ 
+ 0.5Ro{R - RoY sin I cos bO,'^ 
+ V0 sin(/ + 9) cos b, 

— Vr cos{l + 9) cos b, 



Vi = UQsinl — Vq cos / 

+ (i? — -Ro)(-Ro cos / — r cos 6)r2o 

+ {R- Ro)^{Rocosl -rcosb)x (6) 

X O-Sf^Q — ri^o cos 6 

+ vecos{l + 9) + viism{l + 9), 

where N is the number of stars used; j is the current star number; Vr is the line-of-sight 
velocity, Vi = 4.74r/x/cos6 and Vb = 4.74r/Xb are the proper motion velocity components 
in the I and b directions, respectively, with the coefficient 4.74 being the quotient of the 
number of kilometers in an astronomical unit and the number of seconds in a tropical 
year; Vj , Vi' are the measured components of the velocity field (data); wl and wj are the 
weight factors; r is the star's heliocentric distance; the star's proper motion components 
III cos 6 and ni cos 6 and Hb are in mas yr~^ and the line-of-sight velocity Vr is in km s~^; 
Uq,Vq,Wq are the stellar group velocity components relative to the Sun taken with the 
opposite sign (the velocity u is directed toward the Galactic center, v is in the direction 
of Galactic rotation, w is directed to the north Galactic pole), we assume Wq to be 7 km 
s~^, because it is poorly determined without invoking the velocity components Vb, Rq is 



the Galactocentric distance of the Sun; R is the distance from the star to the Galactic 
rotation axis, 

^r^cos^6-2i?orcos6cosZ + i?o. (7) 

Qo is the angular velocity of rotation at the distance i?o; the parameters 0!^ and Qq are, 
respectively, the first and second derivatives of the angular velocity. To take into account 
the influence of the spiral density wave, we used the simplest kinematic model based on 
the linear density wave theory by Lin and Shu (1964), in which the potential perturbation 
is in the form of a travelling wave. Then 

vr = /ijcosx, 
ve = fe sin x, 

where /r and fe are the amplitudes of the radial (directed toward the Galactic center in 
the arm) and azimuthal (directed along the Galactic rotation) velocity perturbations; the 
wave phase x in general form is (Rohlfs 1977) 

X = m[Qp-n(i?)]t + $(i?), (9) 

where for a logarithmic spiral ^{R) — cot(i) In^R/Ro), Eq. (9) for a fixed instant of time 
t can be written as 

X = m[coi{i) HR/Ro) -e]+ xo, (10) 

fn and fg enter into Eqs. (5) and (6) precisely for this approximation; i is the spiral 
pitch angle {i < for winding spirals); m is the number of arms, we take m = 2 in this 
paper; 9 is the star's position angle (measured in the direction of Galactic rotation); Xo 
is the Sun's phase angle, measured here from the center of the Carina-Sagittarius spiral 
arm (R . 7 kpc), as was done by Rohlfs (1977). The parameter A is the distance (along 
the Galactocentric radial direction) between adjacent segments of the spiral arms in the 
solar neighborhood (the wavelength of the spiral density wave). it is calculated from the 
relation 

Am 

tanz = . 11 

2nRo ^ ' 

The weight factors in functional (4) are assigned according to the following expressions 
(for simplification, we omit the index j): 



Wr = So/^S ^ + a^^, ^^^^ 



where 5*0 denotes the dispersion averaged over all observations, which has the meaning 
of a "cosmic" dispersion taken to be 12 km s"^; /3 — ov^jcrvi = 1 is the scale factor that 
we determined using data on open star clusters (Bobylev et al. 2007). The errors of the 
velocities Vi are calculated from the formula 



ay = 4.74r 



\ 




(13) 



The optimization problem (4). (12) is solved for nine unknown parameters Uq^ Vq, Qq, Q'q, 
^0, Jr, fe, i and Xo by the coordinate- wise descent method. 



We estimated the errors of the sought-for parameters through Monte Carlo simulations. 
The errors were estimated by performing 1000 cycles of computations. For this number 
of cycles, the mean values of the solutions essentially coincide with the solutions obtained 
from the input data without any addition of measurement errors. Measurements errors 
were added to such input data as the line-of-sight velocities, proper motions, and distances. 

Here, we take a fixed value of RO based on the review by Foster and Cooper (2010), 
where the weighted mean was Rq = 8.0 ± 0.4 kpc. 

Once the parameters of the rotation curve Q,{R), the perturbation amplitudes and 
fe, and the pitch angle i have been found from the solution of the system of equations 
(5), (6), we are able to determine Qp. For this purpose, we use the approach applied by 
Mishurov et al. (1979). It is based on the assumption about an ellipsoidal distribution of 
residual stellar velocities. According to Lin et al. (1969), 

where A is the amplitude of the spiral density wave potential, — ^1 + is the 

epicyclic frequency, u = m{flp—fl) /nis the frequency with which a test particle encounters 
the passing spiral perturbation, k — m - cot{i)/Ro is the radial wave number, FI)^\x) and 
F^'^\x) are the reduction factors, which are functions of the coordinate x = k'^a\/K^,, 
where ctr is the semimajor axis of the velocity ellipsoid. An expanded form of Eqs. (14) 
and (15) can be found in Mishurov et al. (1979). 

Once Eqs. (14), (15) have been solved, we are able to estimate the ratio of the radial 
component of the gravitational force produced by the spiral arms to the total gravitational 
force of the Galaxy, /^o,, based on the well-known relation (Fernandez et al. 2008; Bobylev 
et al. 2011) 

m 

Our proposed approach to estimating the spiral pattern speed consists in finding the shifts 
in phase Xe from several samples of stars with an age difference At and to determine Jlp 
from relation (9) for the known fl{Ro) 

An^n,- n{R,) = (i7) 

where the phase difference Axq is in radians, and the age difference At is in Myr. Note 
that the change in the Sun's phase does not depend on the Galactic orbit of the Sun; we 
just find all of the parameters being determined for R — Rq and, therefore, it would be 
more appropriate to call Xq the phase of the solar circle. 



2.2 Separate Approach 

In the separate approach, a periodogram analysis based on the Fourier transform to 
determine the periodicities in the velocity components Vr and AVrot is used (Bobylev et 
al. 2008; Bobylev and Bajkova 2010). In this paper, we apply the method of Fourier 



analysis that takes into account the logarithmic pattern of the spiral density wave and 
the distribution of position angles 6 (Bajkova and Bobylev 2012). . 

Initially, we define the heliocentric components of the Cepheid space velocities U and 
V via the observed velocities Vr,Vi, and H based on the well-known relations 

U — cos Z cos 6 — V; sin Z — Vfe cos Z sin 6, 

(18) 

V — Vr sin / cos b + Vi cos I — Vb sin I sin b, 

Subsequently, we find two velocity components: the radial velocity Vr directed from 
the Galactic center to the object and the tangential velocity Vgi in the direction of Galactic 
rotation: 

Ve = Usm9 + {Vo + V)cos9, 

Vr = -C/cos^ + (yo + 1^)sin^, ^ ' 

where Vq = |-Ro^^o| and the position angle 9 is calculated as tan^ = y/{Ro — x), where 
X and y are the heliocentric rectangular coordinates of the stars. Finally, based on the 
velocities Vg, we form the residual tangential velocities of the Cepheids AVrot from which 
the rotation curve found above was subtracted. 

In contrast to the previous method, here we do not assume that the wavelength A is 
the same for the velocity perturbations in the radial and tangential directions and that 
their phases Xq differ exactly by 7r/2. 

Below, we describe the method of spectral analysis of the radial velocity perturbations 
for objects both in the linear approximation of the dependence of the argument in (10) 
on Galactocentric distance and in its exact expression, i.e., by taking into account the 
logarithmic dependence on Galactocentric distance and position angle 9. 

Let there be a series of measured velocities Vr^ (these can be both radial, Vr,, and 
tangential, AVrot, velocities), n = 1, . . . , N, where is the number of objects. The goal 
of the spectral analysis is to separate the periodicity from the data series in accordance 
with model (8)-(ll), which describes a spiral density wave with parameters fR,X{i), and 

Xe- 

Given relation (11) between the spiral pitch angle i and wavelength A for \R—Ro\ ^ Ro 
and small 9, the linear approximation for the logarithm of the argument in (10) can be 
represented as 

In this case, for our harmonic analysis of the velocities, we can apply the standard Fourier 
transform 

_ 1 Af 277- 

^A, = T7 E exp(-j-(i?„ - i?o)), (20) 

n=l '^k 

where V\^ is the kih harmonic of the Fourier transform, Vr^ are the velocity measurements 
for objects with Galactocentric distances Rn,n — 1,2, N, and A^ is the wavelength of 
the kth harmonic, which is equal to D/k, where D is the period of the series being 
analyzed. 

Since we are interested only in the perturbation power spectrum (periodogram) iVx^P 
— , Eq. (20) can be simplified as follows: 

_ I N 27J- 

^A, = T7 E ^Rn exp(-j— 

-'^ n=l 



Note that the derived hnear approximation is acceptable only for analyzing the perturba- 
tions in a small solar neighborhood, while an exact realization of relation (10) is required 
for a wide scatter of distances and position angles for the objects. 

Let us analyze the perturbations as a periodic function of the logarithm of the Galac- 
tocentric distances, for the time being, without allowance for the position angles of the 
objects: 

1 ^ IttR 

^A. = ^ E M-j^ HRn/Ro)). (21) 

-'^ n=l '^fe 

Obviously, if we make the change of variables 

< = ]n{RjR,)R,, (22) 
is reduced to the standard Fourier transform 

K = jfY.V,.M-i^y (23) 

-'^ n=l 

To take the position angles of the objects into account, we will represent Eq. (10) for the 
phase as 

X = Xi- mO, (24) 

where (given (11)) 

^^^^ln(i?/i?„)+xe. (25) 

Substituting (24) into Eq. (10) for the perturbations at the nth point and making standard 
trigonometric transformations, we will obtain: 

Vr„ = /iiCOs(xi„ - men) 

= fn cos xin cos m9n + fn sin xi„ sin m9n (26) 
= fn cos xi„ (cos m9n + tan xi„ sin m6ln) . 

Let us designate 

Vr' = fncosxi, (27) 

It then follows from (26) that 

^Rn = ^< (cos men + tan xi„ sin m^^) . (28) 
Using Eq. (28), let us form a new data series 

^< ^ ^ii„/(cos m9n + tan xi„ sin m9n), (29) 

to which a Fourier analysis can be applied in accordance with (23). A similar relation can 
also be derived for the tangential velocities AVrot- 

Thus, taking into account both the logarithmic pattern of the spiral density wave and 
the position angles of the objects, we obtain the following expression for our spectral 
analysis of the perturbations: 

%. = ^EVH;axp(-,H^). (30) 

-'^ n=l 



The numerical algorithm for realizing (30) consists of the following steps: 

1. The initial series of velocities Vr^ is transformed into the series V^/ in accordance 
with (22). 

2. The power spectrum of the derived sequence Vr,' is calculated based on the Fourier 
transform (23) to obtain an estimate of Xmax that corresponds to the peak of the derived 
power spectrum. 

3. A comb of several Xi{i — 1, . . . , K) with a central Xmax is then specified. 
The following iterations are made for each Aj from the specified comb: 

Step 1. The value of Aj and the initial approximation Xe (for example, equal to zero) 
are substituted into Eq. (25) to calculate Xi for ^ach data reading {n = 1, . . . , N). 

Step 2. Using Eq. (29), the series of velocities Vr>' is transformed into the series V^, . 

■n H„ 

This transformation needs to be regularized to avoid the division by numbers close to 
zero. This is done by assigning a threshold number, say, e, and permission for the division 

is given only when the denominator in Eq. (29) exceeds this number. The best e at which 
the significance of the extracted peak in the spectrum reaches its maximum as a result of 
the iterations at the minimum residual between the solution and the data can be found 
by an exhaustive search for e from some interval. The typical values of e found on model 
problems he within the range [0.01, 0.3]. 

Step 3. The power spectrum of the derived sequence V„, is calculated based on the 

Fourier transform (30) to obtain a new estimate of Xe corresponding to a fixed Aj of the 
derived power spectrum. 

Step 4. The return to the first step is made until the process will converge or diverge. 

Step 5. If the process converged, then we fix the specified Aj and the derived if it 
diverged, then we take the next value Aj+i from the specified comb and make iterations 
1-4 until the value of A at which the process converges will be found. 

4. The power spectrum is calculated for the values of A and x& found based on Eq. 
(30) with the goal of a further analysis. 

A more detailed description of the algorithm, the results of its testing on model data, 
and the questions of estimating the significance of the extraction of peaks in the peri- 
odogram and estimating the errors in the spiral density wave parameters are given in 
Bajkova and Bobylev (2012). 



3 RESULTS AND DISCUSSION 

3.1 Results of the Simultaneous Solution 

The results of the solutions of the system of equations (5) and (6) obtained from four 
samples of Cepheids are presented in Table 2: for Cepheids with periods longer than 9 
days, i.e., the youngest ones in our sample, in the first column; for Cepheids with periods 
from 5 to 9 days, i.e., middle-aged ones, in the second column; for Cepheids with periods 
shorter than 5 days, i.e., the oldest ones, in the third column; the solution obtained from 
all Cepheids is given in the fourth column. 

The parameters of the Galactic rotation curve Qq, Qq, and Qq, calculated from the 
sample of middle-aged Cepheids are in good agreement with the results of analyzing blue 
supergiants (Zabolotskikh et al. 2002), OB associations (Mel'nik and Dambis 2009), 
and Galactic masers (Bobylev and Bajkova 2010; Stepanishchev and Bobylev 2011). In 



Table 2: Kinematic parameters found from Cepheids 



Parameters 


7 ~) \ f\(i 

> 9 


rd ^ I ~> r\(i 

5 < -T < 9 


F < b 


A 11 

All 


Uq, km s 








0.9 ± 1.2 


7.5 ± 0.9 


8.0 ± 1.3 


7.0 ± 0.5 


Vq, km s 


-1 






11.6 ± 0.8 


10.9 ± 0.6 


15.6 ± 1.3 


11.7 ± 0.3 


\Lq, km s 


_i 


kpc 


-I 


26.1 ± 0.9 


30.4 ± 1.0 


24.1 ± 1.2 


27.5 ± 0.5 


km s 




kpc'^ 




—3.95 ± 0.13 


—4.34 ± 0.13 


—4.72 ± 0.25 


—4.12 ± 0.10 


iZg, km s 


-1 


kpc 


-3 


0.79 ± 0.10 


0.69 ± 0.14 


2.44 ± 0.28 


0.85 ± 0.07 


/«, km s 








—9.8 di 1.6 


O C _l_ 1 1 

—8.5 ± 1.1 


IOC _l_ 1 T 

— iz.D ± i.7 


G. o _i_ r\ T 
—0.8 di (J. 7 


Je, km s 


1 






1 O 1 1 A 

1.3 ± 1.9 


2.7 ± 1.1 


8.2 ± 1.4 


3.3 ± 0.5 


deg 








—5.2 ± 0.3 


—4.0 ± 0.1 


—6.6 ± 0.5 


—4.6 ± 0.1 


X0, deg 








— 14s ± 14 


1 no n 
—193 ± 9 


—234 ± 10 


—193 ± 5 


(7o, km s~ 


-1 






1 o o 
Iz.z 


1 o 
Iz.o 


Iz.o 


1 Q /I 

lo.4 


A, kpc 








2.3 ±0.5 


1.8 ±0.1 


2.9 ±0.4 


2.0 ±0.1 










61 


72 


52 


185 


h, Myr 








33.9 


59.2 


84.9 




^2, Myr 








55.4 


95.2 


135.0 




o-R, km s- 


-1 












14 


r2p, km s" 


-1 


kpc" 


-1 








23.5 


frO 














0.04 ± 0.01 



Note: is the number of Cepheids in the sample, cxo is the error per unit weight obtained 
when solving the system of equations (5), (6), or is the semimajor axis of the velocity ellipsoid, 
Op was estimated by the method of Mishurov et al. (1979). 

these papers, the angular velocity of Galactic rotation (^o is ~30 km s~^ kpc~^. For 
unknown reasons, the youngest Cepheids revolve around the Galactic center with a lower 
velocity. Note that and l^g, that we found based on the entire sample of Cepheids 
(the last column in Table 2) are in good agreement with the results of the analysis of 
only the proper motions of Cepheids from the Hipparcos catalog performed by Feast and 
Whitelock (1997), where = 27.2 ± 0.9km s"^ kpc'^ was found. 

It can be seen from Table 2 that the amplitudes of the velocity perturbations in the 
radial direction, f^,, differ significantly from zero in all cases, while in the tangential 
direction, Jq,, they are significant for the sample of old Cepheids (and for the entire 
sample) . The relationship between the amplitudes found is in agreement with the results of 
the analysis of blue supergiants performed by Zabolotskikh et al. (2002), = —6.6 ±2.5 
km s~^ and fg = 0.4 ± 2.3 km s^^ for m = 2, and a sample of young Cepheids with 
similar values of these parameters. Previously (Bobylev and Bajkova 2011), we found 
fR = -12.5 ± 1.1 km s-^ fe = 2.0 ± 1.6 km s-^ and i = -5.3 ± 0.3° for m = 2 with 
Xe — ~91 i 4° from data on 0B3 stars with an independent distance scale determined 
from interstellar Ca II absorption lines. 

In principle, the Sun's phases found (Table 2) are consistent with the results of an- 
alyzing various samples of Cepheids, for example (with the phase measured from the 




Figure 2: Radial velocities for the sample of young Cepheids (a) and their tangential 
velocities (b) versus Galactocentric distance R; the vertical dotted lines mark the Sun's 
position, the solid curve in panel (b) indicates the rotation curve, the lines with long 
dashes mark the la confidence intervals. 



Carina- Sagittarius arm): Xq = —165 ± 1° (Byl and Ovendcn 1978), Xq = —150° from 
red supergiants and Cepheids (Mishurov et al. 1979), Xq = —290 ± 16° (Mishurov et al. 
1997), and Xq = -320 ±9° (Mishurov and Zenina 1999) from relatively old Cepheids with 
periods P < 9'^. 

In Figs. 2 and 3, the radial, Vr, and tangential, Vrot = Ve, velocities (calculated 
from Eqs. (19)) are plotted against the Galactocentric distance R for two samples of 
Cepheids. young and middle-aged ones. In each case, the rotation curve was constructed 
in accordance with the results of Table 2. A wave structure in the radial velocities Vr is 
clearly seen in Figs. 2 and 3. 

As can be seen from Table 2, the phase Xe changes noticeably with the mean age of 
the sample. We also see that the mean Cepheid ages calculated using calibrations (2) 
and (3) differ by a factor of 1.5. Applying Eq. (17) gives |Af2|(i_2) = 9.9 km s"^ kpc~^ 
(columns 1, 2, 3 in Table 2), |A1]|(2_3) = 9.4 km s"^ kpc^^ and |Al]|(i_3) = 9.0 km s'^ 
kpc~^ for the calibration of Efremov (2003) (Eq. (3)). In this case, the mean value of the 
difference |Af2| = 9.4 km s~^ kpc~^, then Qp = 18.1 km s~^ kpc~^ (for Qq = 27.5 ±0.5 km 
s^^ kpc^^ found from the entire sample of Cepheids), which is in satisfactory agreement 
with the result obtained from Eqs. (14) and (15) by the method of Mishurov et al. (1979), 
Qp — 23.5 km s~^ kpc~^. 

The analogous values found from the age cahbration of Bono et al. (2005) (Eq. (2)) 
are considerably larger: |Ar2|(i_2) = 15.4 km s~^ kpc~^, |Af2|(2-3) = 14.7 km kpc~^ 
and |Ai7|(i_3) = 14.0 km s'^^ kpc~^. In this case, the mean value of the difference \AQ\ = 
14.7 km s~^ kpc~^, then Qp = 12.8 km s~^ kpc~^, which is in much poorer agreement 
with other known data. 




Figure 3: Radial velocities for the sample of middle-aged Cepheids (a) and their tangential 
velocities (b) versus Galactocentric distance R. The designations are the same as those in 
Fig. 2. 



Comparison of the results obtained allows us to opt for the age calibration of Efremov 
(2003). 

The estimate of fro = 0.04 ±0.01 that we obtained from Eq. (16) is much smaller than 
fro = 0.15, found from a sample of 192 long-period Cepheids by Mishurov et al. (1979). 
At the same time, our estimate is in good agreement with f^o — 0.05, suggested by Yuan 
(1969). 

3.2 Results of the Separate Approach 

Figure 4 presents the Galactocentric radial velocities Vr of the Cepheids and the corre- 
sponding power spectra. The distance R indicated in the plots was calculated from the 
relation R = \n{R/Ro)Ro. which takes into account the logarithmic pattern of the spiral 
density wave (22). In contrast to Figs. 2 and. 3, the data here are presented in the form 
of "pulses" without indicating their random errors. 

As in the previous case (Table 2), this approach failed to reveal a significant wave 
in the tangential velocities for the sample of the youngest Cepheids. The proposition of 
the density wave theory that the perturbations in the radial direction propagate faster 
than those in the tangential one can serve as an explanation of this fact. Therefore, the 
youngest stars have not yet had time to respond to the perturbations in the tangential 
direction, while the situation for the older stars gradually levels off. 

Figure 5 shows the residual tangential velocities AVrot for the sample of middle-aged 
Cepheids and their power spectrum. 

We clearly see from Figs. 4 and 5 and the data of Table 2 that the wavelength A 




Figure 4: Radial velocities for the sample of young Cepheids (a) and their power spectrum 
(b), middle-aged Cepheids (c) and their power spectrum (d), and old Cepheids (e) and 
their power spectrum (f); the vertical dotted lines mark the Sun's position. 




Figure 5: Residual tangential velocities for the sample of middle-aged Cepheids (a) and 
their power spectrum (b); the corresponding radial velocities are shown in panels (c and 
d) of Fig. 4. 



changes, depending on the sample and the velocity character. For example, the difference 
in A calculated from the radial and tangential velocities of middle-aged Cepheids is about 
0.4 kpc. Comparison of the data in Fig. 5 with panels (c) and (d) in Fig. 4 leads us 
to conchidc that there is a phase shift close to 7r/2 (the exact value is 110°), which is in 
agreement with the prediction of the linear density wave model. It can be seen from the 
sample of middle-ages Cepheids that the Sun is equidistant from the centers of the two 
nearest segments of the spiral arms. 

For a more reliable determination of the mean AQ, we invoked our previously obtained 
data on the kinematics of young 0B3 stars (Bobylev and Bajkova 2011), where the Sun's 
phase Xq = — 91 ± 4° was obtained. For these stars, we took a mean age of 8 Myr. 
We calculated the mean \ AQ\ = 9.8 ± 0.6 km s~^ kpc~^ using six independent differences 
formed from four values of Xe- —91°, —148°, —193° and —234°. We used the mean Cepheid 
ages calculated from Efremov's cahbration. Here, the error {eAn)stat — 0.6 km s""^ kpc~^ 
was calculated from the convergence of the results. We can also estimate the systematic 
error of the method, {sAn)syst, from the relation 

for the adopted typical values: = 10° ■ (7r/180°), = 50° ■ (7r/180°), At = 50 Myr 
and Et = 10 Myr, then {sAnjsyst = 2.8 km kpc~^. As a result, we have \Afl\ = 
9.8 ± O.Gstat ± ^.Ssyst km s'^ kpc'^ 

Note that we found x& — — 130 ± 10° from a sample of masers in regions of active star 
formation (Bobylev and Bajkova 2010). Thus, they are intermediate between the 0B3 
stars and young Cepheids. Although there are no estimates of their individual ages, there 
are very massive O stars as well as high- and low-mass proto stars among them. The 



paradoxical value of x& = —130° (the expected Xe ~ —90°) can be explained by the fact 
that the kinematics of these recently formed stars probably reflects considerably earlier 
stages in the Galactic motion of the regions of active star formation. 

The Galactic spiral pattern speed calculated by our proposed method from the entire 
sample of Cepheids {Qq = 27.5 km s"^ kpc~^) is Qp = 17.7 km s~^ kpc~^. If, however, 
we take Qo = 30 km s~^ kpc~^ known from our analysis of masers and 0B3 stars, then 
we obtain Qp — 20.2 km s~^ kpc~^. According to these data, the co- rotation circle in 
the Galaxy is located at a distance from R —10 to 12 kpc. These estimates are valid 
for a two-armed Galactic spiral pattern (m = 2). In the opinion of several authors, the 
Galactic spiral pattern is a four-armed one (m = 4). In this case (see Eq. (17)), the 
difference we found will be — r2o| ~ 5 km s^^ kpc~^. Then, the spiral pattern speed 
will be Qp — 22.5 km s'^ kpc"^ from the entire sample of Cepheids (fio ~ 27.5 km s~^ 
kpc~^) and flp — 25 km s""^ kpc~^ for Qq — 30 km s"^ kpc""^; accordingly, the co-rotation 
circle will be still closer to the Sun. 

4 CONCLUSIONS 

We analyzed the space velocities of Galactic Cepheids to study the peculiarities of the 
Galactic spiral density wave. For this purpose, we used 185 Cepheids with proper motions 
mainly from the Hipparcos catalog and line-of-sight velocities from various sources. We 
divided the entire sample into three parts, depending on the pulsation period, which 
reflects well the mean Cepheid age. 

First, based on the entire sample of Cepheids and taking Rq = 8 kpc, we found the 
components of the peculiar solar velocity {uq,Vq) = (7.6, 11.7) ± (0.5,0.3) km s~^, the 
angular velocity of Galactic rotation = 27.5 ± 0.5 km s~^ kpc~^ and its derivatives 
Q'q = -4.12 ± 0.10, km s"^ kpc'^ and ^Iq = 0.85 ± 0.07, km s"^ kpc"^ the amplitudes of 
the spiral density wave = —6.8 ±0.7 km s^^ and fe = 3.3 ±0.5 km s^^, the pitch angle 
of a two-armed spiral pattern i = —4.6 ± 0.1° (then A = 2.0 ±0.1 kpc), and the phase of 
the Sun in the spiral density wave x& = —193 ± 5°. 

The Cepheids with pulsation periods P from 5 to 9 days (middle-aged) (Qq = 30.4±1.0 
km s~^ kpc~^ for them) show the fastest Galactic rotation. The amplitude of the radial 
velocity perturbations fji caused by the spiral density wave in each of the three samples is 
about 9 km s^^; the amplitude of the tangential velocity perturbations fg increases from 
zero for the youngest Cepheids to ~4 km s^^ for older ones. 

We found that the Sun's phase x& changes noticeably with the mean age of the sample. 
For a more accurate estimation of the phase, we applied a Fourier analysis of the Cepheid 
radial velocities, while for middle-aged Cepheids we also managed to reliably determine 
the parameters of the tangential velocity perturbations by this method. 

From our analysis of the phase shifts, we determined the mean value of the angular 
velocity difference AQ = — Q, which depends significantly on the calibrations used to 
estimate the individual ages of Cepheids. For a more reliable determination of the mean 
AQ, we invoked our previously obtained data on the kinematics of young 0B3 stars. As 
a result, we calculated the mean AQ using six independent differences. 

We showed that when the individual ages of Cepheids derived from the period-age 
calibration of Bono et al. (2005) is used, the differences Qp — Qq exceed those derived 
from the cahbration of Efremov (2003) by a factor of 1.5-2. Simultaneously, we obtained 



an estimate of Qp = 23.5 km kpc~^ from our analysis of the reduction factors by 
the method of Mishurov et al. (1979). As a result, this allowed us to opt for Efremov's 
calibration , using which we found |Qp — QqI — 9.8 ± O.dgtat i '^■^syst km s~^ kpc"-*^. 

The Galactic spiral pattern speed that we calculated by our proposed method based 
on the entire sample of Cepheids {Qq = 27.5 km s~^ kpc~^) is Qp = 17.7 km s~^ kpc~^. If, 
however, we take flo = 30 km s~^ kpc~^ known from our analysis of masers and 0B3 stars, 
then we will obtain Qp — 20.2 km s~^ kpc~^. According to these data, the co-rotation 
circle in the Galaxy is located at a distance from R = 10 to 12 kpc. 

Based on the entire sample of Cepheids, we estimated the ratio of the radial component 
of the gravitational force produced by the spiral arms to the total gravitational force of 
the Galaxy to be fro = 0.04 ± 0.01. 

5 ACKNOWLEDGMENTS 

We are grateful to the referees for valuable remarks that contributed to a significant 
improvement of the paper. This work was supported in part by the "Nonstationary 
Phenomena in Objects of the Universe" Program of the Presidium of the Russian Academy 
of Sciences and grant no. NSh-1625.2012.2 from the President of the Russian Federation. 
In our work, we used the SIMBAD search database and the DDO Cepheid database. 



6 REFERENCES 

1. LA. Acharova, Yu. N. Mishurov, and V. V. Kovtyukh, Mon. Not. R. Astron. Soc. 

420, 1590 (2012). 

2. A. T. Bajkova and V. V. Bobylev, Astron. Lett. 38, 549 (2012). 

3. L. N. Berdnikov, A. K. Dambis, and O. V. Vozyakova, Astron. Astrophys. Suppl. Ser. 
143, 211 (2000). 

4. V. V. Bobylev, A. T. Bajkova, S. V. Lebedeva, Astron. Lett. 33, 720 (2007). 

5. V. V. Bobylev, A. T. Bajkova, and A. S. Stepanishchev, Astron. Lett. 34, 515 (2008). 

6. V.V.Bobylev, andA.T. Bajkova,Mon. Not.R. Astron. Soc. 408, 1788 (2010). 

7. V. V. Bobylev and A. T. Bajkova, Astron. Lett. 37, 526 (2011). 

8. V. V. Bobylev, A. T. Bajkova, A. Myllari, and M. Valtonen, Astron. Lett. 37, 550 
(2011). 

9. G. Bono, M. Marconi, S. Cassisi, et al, Astrophys. J. 621, 966 (2005). 

10. J. Byl and M. W. Ovenden, Astrophys. J. 225, 496 (1978). 

11. Yu. N. Efremov, Astron.Rep. 47, 1000 (2003). 

12. M. Feast and P.Whitelock, Mon. Not. R. Astron. Soc. 291, 683 (1997). 

13. D. Fernandez, F. Figueras, and J. Torra, Astron. Astrophys. 480, 735 (2008). 

14. T. Foster and B. Cooper, ASP Conf. Ser. 438, 16 (2010). 

15. P. Fouqu, P. Arriagada, J. Storm, et al., Astron. Astrophys. 476, 73 (2007). 

16. O. Gerhard, Mem. Soc. Astron. Ital. Suppl. 18, 185 (2011). 

17. G. A. Gontcharov, Astron. Lett. 32, 759 (2006). 

18. E. Hog, C. Fabricius, V. V. Makarov, et al, Astron. Astrophys. 355, L27 (2000). 

19. F. van Lccuwcn, Astron. Astrophys. 474, 653 (2007). 

20. C. C. Lin and F. H. Shu, Astrophys. J. 140, 646 (1964). 

21. C. C. Lin, C. Yuan, and F. H. Shu, Astrophys. J. 155, 721 (1969). 



22. A. M. Mel'iiik, A. K. Dambis, and A. S. Rastorguev, Astron. Lett. 25, 518 (1999). 

23. A. M. MePnik and A. K. Dambis, Mon. Not. R. Astron. Soc. 400, 518 (2009). 

24. Yu. N. Mishurov, I. A. Zenina, A. K. Dambis, et al, Astron. Astrophys. 323, 775 
(1997). 

25. Yu. N. Mishurov, E. D. Pavlovskaya, and A. A. Suchkov, Sov. Astron. 23, 147 (1979). 

26. Yu. N. Mishurov and I. A. Zenina, Astron. Astrophys. 341, 81 (1999). 

27. M. E. Popova, Astron. Lett. 32, 244 (2006). 

28. A.S.Rastorguev,E. V. Glushkova, A. K. Dambis, and M. V. Zabolotskikh, Astron. 

Lett. 25, 595 (1999). 

29. K. Rohlfs, Lectures on Density Wave Theory (Springer, Berhn, 1977). 

30. A. S. Stepanishchev and V. V. Bobylev, Astron. Lett. 37, 254 (2011). 

31. C. Yuan, Astrophys. J. 158, 889 (1969). 

32. M. V. Zabolotskikh, A. S. Rastorguev, and A. K. Dambis, Astron. Lett. 28, 454 
(2002). 

33. N. Zacharias, C. Finch, T. Girard, et al., CDS Strasbourg, 1/315 (2009) 

34. The Hipparcos and Tycho Catalogues, ESA SP-1200 (1997). 



